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A Langevin equation for the complex amplitude of a single-mode Bose-Einstein condensate is 
derived. The equation is first formulated phenomenologically, defining three transport parameters. 
It is then also derived microscopically. Expressions for the transport parameters in the form of 
Green-Kubo formulas are thereby derived and evaluated for simple trap geometries, a cubic box 
with cyclic boundary conditions and an isotropic parabolic trap. The number fluctuations in the 
condensate, their correlation time r c , and the temperature-dependent collapse-time of the order 
parameter as well as its phase-diffusion coefficient are calculated, 



I. INTRODUCTION 

Bose-Einstein condensation in a weakly interacting Bose-gas in three dimensions in the thermodynamic limit of an 
infinitely extended system is a second order phase transition in which an order parameter, the macroscopic wave- 
function, appears spontaneously with a fixed but arbitrary phase, turning the global £/(l)-gauge-symmetry connected 
with particle-number conservation into a spontaneously broken or hidden symmetry. The rigidity of the phase of the 
order parameter against local perturbations and the absence of any phase diffusion gives rise to the Goldstone modes, 
which take the form of collision-less (zero)-sound or hydrodynamic sound, respectively, depending on whether the 
sound frequency is in the collision-less mean-field regime or in the collision-dominated regime Q . 

In finite systems, and thus also in all trapped Bose-gases, sharp phase-transitions are impossible and hidden sym- 
metries in a rigorous sense cannot appear. However, a macroscopic wave-function describing a Bose-Einstein con- 
densate still exists ||. Its phase cannot be stable and must undergo a diffusion process, which restores the U(l) 
gauge-symmetry over sufficiently long time intervals. This diffusion process is different from the Goldstone modes 
mentioned before, which are oscillations around a fixed value of the phase and do not restore the symmetry. Rather 
the Goldstone modes show up either as collision-dominated hydrodynamic phonons or as collision- less phonons, which 
have also been observed in the finite Bose-Einstein condensates. In the present paper I would like to discuss the 
dynamics of the complex amplitude of a Bose-Einstein condensate containing a finite number of particles, and in par- 
ticular analyse the diffusion of its phase. My discussion will extend and correct in several respects the work published 
in | . 

The stability of the phase-difference between the macroscopic wave functions of two Bose-Einstein condensates 
in a trap has been measured. In the experimental set-up the relative phase was measured using a time-domain 
separated-oscillatory-field condensate interferometer. Over the time-interval of 100 ms scanned in the experiment the 
relative phase was found to be robust. This experimental result demonstrates that the macroscopic wave functions 
of the condensates cannot be considered as quantum mechanical wave-functions of many particle systems entangled 
with each other, whose decoherence would indeed be extremely rapid. Rather, the macroscopic wave functions 
are appropriately viewed as robust classical objects, their quantum mechanical origin (just like magnets, crystals 
etc) notwithstanding. This does, of course not preclude that there may be quantum effects, for finite condensates, 
which lead to corrections of the dynamics described by the underlying classical wave-equation, the well-known Gross- 
Pitaevskii equation |tJ . In a number of papers || the dispersion of the phase of a trapped Bose-Einstein condensate 
at zero-temperature was considered, which is due to fluctuations Sfi of the chemical potential [i in a finite system with 
fixed particle number. An extension of this mechanism to finite temperature has also been proposed |^]. This effect is 
not an irreversible phase diffusion but corresponds to an effect of inhomogeneous broadening, similar to the de-phasing 
of processing spins occurring in spin-systems due to inhomogeneous broadening. Like the decay of the magnetization 
can be reversed in spin-echos, the decay of the order parameter expectation value in Bose-Einstein condensates due 
to a finite variance of Sfi is in principle reversible in 'revivals'. Experiments in Bose-Einstein condensation are done 
at temperature fcgT 3> Tiui and often even at ksT 3> /i, where uj is the geometrical mean of the three main trap 
frequencies. A process of phase-diffusion process should occur in such a regime due to the interaction of the condensate 
with a thermal bath of collective modes and quasiparticles. An estimate of this phase-diffusion is of interest for the 
theory of atom lasers, because the fundamental limit of the line-width of an atom laser for a given temperature 
depends on it similar to the 'Schawlow-Townes'-formula jl(| for the line-width of a laser. 

In this paper a theory of dissipation and thermal fluctuations of a trapped Bose-Einstein condensate will be formu- 
lated. First a phcnomcnological framework for the theory in the form of a Langevin equation will be given in which 
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dissipation appears via a phenomenological parameter and the fluctuation-dissipation relation is invoked to relate it to 
maximally three intensity-coefficients of the fluctuations. The solution of the Langevin-equation then determines the 
relaxation of the condensate-number and the diffusion of the phase, quite similar to the dynamics of a laser-amplitude 
above threshold. Then the Langevin equation is derived from the microscopic theory and formulas for the phenomeno- 
logical parameters are derived. These are evaluated for a box-like trap and an isotropic harmonic trap-potential as a 
function of temperature, particle-number and scattering length. The final section contains a discussion of our results 
and a comparison with earlier related work. The theory presented here makes no claim to apply to the critical regime, 
nor can we examine here to what extent it covers the regime below but close to T c , where it may be important to 
take the dynamics of the thermal cloud of non-condensed atoms into account in addition to the excitations from the 
condensate. 



II. MICROSCOPIC EQUATIONS OF MOTION 

The weakly interacting Bose-gas in a trap in standard notation is described by the Hamiltonian 

i xi>+{-^V 2 + V(x)-{ t x) + ^+i>}i>. (2.1) 

The total number of atoms N is fixed, i.e. the Hilbert space is the restriction of the Fock-space of ip to the subspace 
on which N = N is satisfied, (fj) is the average of the chemical potential, which is a fluctuating quantity in a system 
where N is fixed. Later-on we shall denote the fluctuating part of the chemical potential with A/x. The presence of a 
Bose-Einstein condensate in equilibrium means that many (No 3> 1) particles occupy a single mode of a macros copi c 
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classical matter wave, determined as the mode of lowest energy of the classical Hamiltonian corresponding to eq. 
The latter is obtained by replacing in H the field operator ip(x) by the classical field i/j(x) = y/ (No) exp(i<p)ipo(x). We 
shall restrict our attention to sufficiently low temperatures below the critical temperature T c so that the interaction 
of the condensate with the mean field of the thermal cloud of non-condensed particles is negligible, fn this way one 
finds that the condensate mode ipo( x ), which we take to be normalized to f, satisfies the Gross-Pitaevskii equation 

§ 

- (h 2 /2m)V 2 ^o + (V(x) + Uo(N )\i>o(x)\ 2 )i>o = Wo- (2.2) 
For given (No) the average value of the chemical potential /i follows by imposing the normalization condition 

d 3 x\$o(x)\ 2 = 1 (2.3) 

on the solution of the Gross-Pitaevskii equation, and thereby (fx) like ipo(x) becomes a function of the mean atom 
number in the condensate (No). As an important consequence of this fact the chemical potential of the system can 
be expressed as a function of the average number of atoms in the condensate alone. (No) differs from N, the fixed 
total number of atoms, by the average number (A 7 ) of non-condensed atoms, which needs to be calculated for given 
(No). The condition N — (No) + (N 1 ) then fixes (No) self-consistently. In the experimentally realized Bose-Einstein 
condensates it is possible to measure (No) directly with reasonable accuracy as a function of temperature, and in 
practice it is therefore reasonable to regard (No) as an experimentally given and known function of temperature. The 
space-dependent mean number density of the condensate is no(x) = (No)\ipo(x)\ 2 . We shall take the mode function 
ipo(x) in the Gross-Pitaevskii equation as real and positive. (This also means we are not considering condensates 
containing vortices) . The physical phase of the condensate is not carried by its mode-function ipo but by its complex 
amplitude denoted as ceo, where ao = y/No expicp. 

If |ao| 2 makes a small fluctuation away from its equilibrium value (No) the condensate mode function ipo will no 



longer satisfy eq.(2.2) but will change its form slightly. We shall assume that such fluctuations of A — \ceo\ 2 occur 
on a sufficiently large time-scale, that the new form is again determined by the Gross-Pitaevskii equation, but for the 
changed condensate number |«o| 2 and a correspondingly changed chemical potential /io determined uniquely by |ao| 2 , 
i.e. in eq.(2.2) the replacements (ipo, (No), ((J,)) — > (ipo, \ao\ 2 ,Ho) have to be made in this case, 



- (H 2 /2m)V 2 iPo + (V(x) + Uo\a \ 2 \M x )\ 2 )^o = MoV>o- (2.4) 

We cannot expect, in general, that in any given nonequilibrium state the difference defined by Ao/i = /-to — (n) is 
the total deviation of the chemical potential from its equilibrium value, because there may obviously be states with 
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|cko| = (No) which differ in other respects from the equilibrium state and may therefore have /i 7^ (fx). Therefore we 
use the notation hq for the part of the nonequilibrium chemical potential determined by |cco| 2 - 

The presence of the highly occupied condensate mode makes the decomposition of the Heisenberg field-operator 

tj)(x, t) = (\a \ exp (i<l>)ip (x) + x(x, t j) cxp (-i(/j,)t/n) (2.5) 

useful, where we follow Bogoliubov and describe the condensate classically. x(x, t) is taken to be the field operator 
for the particles outside the condensate. We shall assume that the temporal changes in (j> can be considered as slow 
on the time-scales of the dynamics of X- The phase (j> an d amplitude |ao| are additional c- number variables in ( |2.5| ). 
Therefore the taking of expectation values has from now on to include an integration over a distribution of |ao| and 
in addition an integration over all values of (j>. Since the total number N is fixed (ip) = must hold for all times. 
However, it will also be useful to consider expectation values in the Fock-space of the operators x> X + alone without 
averaging over <\>. Such expectation values will be denoted as (...)</,. 

Gauge invariance, strictly speaking, is lost by splitting off a c-number term from the field-operator. However, this 
symmetry is saved by adopting the rule that the phase <j) of the c-number term in the decomposition also changes 
under a gauge transformation according to <p — > <p + e. By this device we take into account the fact that the same 
change of phase would have occurred automatically, if we had not replaced the condensate term by a c-number. The 
generator of gauge transformations is thus taken as 



d 

N = i — 



+ / d A x X +X, (2.6) 



from which it is clear (cf. eq.( 2.24 )) that i-§^\ x ^ x + ^ s a representation of A^o-Q The canonical-conjugate is the phase cj) 
with 

exp(i0) = exp(d/dN )£,x+ ■ (2.7) 

Via (^^) the Hamiltonian furthermore splits up according to H = Hq + H\ + H2 + H3 + H4 where H n comprises 
the terms of H which are of n-th order in x, x + . Explicitly 



//(. : : ^o| 2 yd 3 ^o{-i:V 2 + y(a ; )- / i + ^|ao| 2 |^o| 2 Wo + (Mo-(M»|ao| 2 (2.8) 



//, : |n„| / d 3 a:{e-^x(-|-V 2 + V(x) - (//) + {/oKlVo^o + (&•*)} (2.9) 

2m 



r h 2 

H 2 = J d?x{x + {-—V 2 + V{x)~^)x 



El 
2 

+ (Mo - (v))x + x} 



+ f\ a ^o( e ~ 2 ^ 2 + ^X +2 + 4X+X) (2.10) 



U \a \ f d 3 x^oX + (e~^X + ^X + )x (2.11) 



H i= El I d 3 x X + X + XX- (2.12) 



Ha 
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Using the Gross-Pitaevskii equation (|2.4|) and its derivative with respect to |ao| 2 we derive in the appendix 



1 This operator with fixed X,X + has to be well distinguished from the unrestricted derivative operator i{d/dcj>), which is a 
representation of the total particle number N and has as formal canonical-conjugate <j) with exp(i^) = exp(<9/cW). 
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hoi 2 

H = I dN (v (N ) - 

'<JV > 



The term Hi can be simplified using the Gross-Pitaevskii equation ( |2.4[ ) and becomes then 

Hi = |a |(Mo - (M» / ^(e^X + e^x + )V>o 



(2.13) 



(2.14) 



This expression will be seen to vanish below due to an orthogonality condition. 

The first part of Hi in (2.10) is diagonalized by introducing quasi-particle operators a„, d+ defined by the standard 
Bogoliubov transformation, with time-dependent <p(i) 



X(x) = ^2 (u u (x)a v + v*(x)a+ 



(2.15) 



The u Vl v v satisfy the usual Bogoliubov- Fetter equations 

( -^J 2 + U eS [x)-Ku v K(x) 

\ K(x) -^v 2 + u eS (x) + hu v 



u v {x) 
v v (x) 



= 0, 



(2.16) 



with the abbreviations 



U eS (x) = V(x) ~ mo + 2U \a \ 2 Mx) 2 
K(x) = \a \ 2 UoMx) 2 - 



The Hamiltonian H2 now takes the form 



H 2 = ^Tioj v {a+a v + \v v \ 2 ) + (fi - (fj) + %4>) J d 3 x X + x 



(2.17) 



(2.18) 



The coefficients u u ,v l/ and the mode frequencies lu u become also functions of |ao| and fluctuate (slowly) with that 
number. Their equilibrium values will be denoted by U v ,v v and Q u and the corresponding operator \ according to 
(2.15) as x. 

Eq. (2.16) is consistent with the ort ho- normality conditions 



d 3 r(u* u v^ - u^Vy) = 



(2.19) 
(2.20) 



which guarantee the Bose commutation relations of the a v , a+. A formal solution of eq. (2.16) at zero energy Tiui v = 
is given by the condensate 



u v {x) = -v* v {x) = tp (x) , uj v = . 



(2.21) 



This solution is obviously not normalizable in the required sense (2.1E) to furnish an acceptable solution for the u„, v v 
and must therefore be excluded from the sum over the terms containing the operators d„,dj . The existence of this 
formal solution implies however that the properly normalizable solutions u Vl v u and the condensate mode ipo satisfy 
the important orthogonality relation 



It follows with ( 2.15 ) that 



d 3 xip (uv +v u ) = 0. 



(2.22) 



(2.23) 



which in turn implies that the reduced expression (2.14) for H± vanishes. Using the property (2.23) one can verify 
that the decomposition (|2.5|) of tfj implies 
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N= |a | 2 + TV' 



(2.24) 



with 



N> = / d 3 xx+(x)x(x) 



(2.25) 



which serves as a definition of N$ = |cto | 2 ■ The mean thermal density n' in equilibrium can now be determined via 



n'(x) = (* + (*)£(*)) = J2{(\M*)\ 2 + \v v {x)\ 2 )n v + \M*)\ 2 } 



(2.26) 



with n v — (exp(/3%u u ) — 1) . 

The fluctuations of Nq are similarly fixed by 



(AiV 2 ) = (iV 2 ) - (N ) 2 = (AN' 2 ) 



^]T<U,(rV + l) f 



d 3 cc(u*(x)w^'(a;) + u„/ 



+{n v n v > +-(n„ + rv + f )) 



d 3 x(u u (x)v U '(x) + u u >(x)v v (x)) 



(2.27) 



(2.28) 



They have been evaluated in JL4| and are also needed below (see eq.( p.20 )). For work in the mathematical physics 
literature on number-fluctuations in the condensate of the ideal Bose-gas and models of the interacting Bose-gas see 
jr2| , fL3[ and references given there. For an alternative proposal to define and calculate the number-fluctuations in a 
Bose-condensate see [[15| . 

After the transformation (2.15) the Hamiltonian is now in the form 

H = H + H 2 +H a + H 4 (2.29) 

with H ,H 2 ,H 3 , Hi given by eqs. (|2~l^pl8|,plpl^) . 



III. L ANGEVIN-EQUATION OF THE CONDENSATE AMPLITUDE 



Neither the Gross-Pitaevskii equation nor the Bogoliubov-Fetter equations furnish an equation for the condensate 
amplitude ao = \/ Nq expiifi. To find such an equation phenomenologically we first turn to a macroscopic quantity like 
the entropy S(\ao\ , N) for a fixed particle number N, but restricted to a fixed arbitrary value of ao — \/No exp («</>), 
where Nq is the instantaneous number of particles in the condensate and different from the equilibrium value (Nq) 
corresponding to the maximum of S(\ao\ 2 ,N). Thus (No) is a function of N. The fluctuations of No in the closed 
system formed by the trapped condensate after the evaporative cooling has been switched off are governed by a 
canonical Boltzmann-Einstein distribution 

P(No) = n- 1 e^p(S(\ao\ 2 ,N)/k B ). 

We shall restrict ourselves to temperatures in the condensed regime outside the critical regime, where (Nq(N)) is 

much larger than its root mean square y/(AN'^(N)) = y(AN> ) = ((N' ) - (N') 2 ) 1 / 2 , which is also a function of N. 
Then S(\ao\ 2 , N), expanded to lowest order around its maximum, takes the form 

S(\a \ 2 ,N) = S^(N) + AS(\a \ 2 ,N) 

with 
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a qn ,2 An , (|«o| 2 - (Np)) 2 . . 

AS(\a \ ,N) = -k B 2 (AN 2 ) ^ ' 

The entropy S'daol 2 , not only determines the equilibrium-distribution of the condensate amplitude, but appears 
also in its equation of motion, both in the conservative part of the dynamics as a conserved quantity, and in the 
dissipative part as a potential for the irreversible part of the dynamics. Let us first consider both parts separately. 

The conservative part of the dynamics of ao is connected with the dynamics of its phase <p. According to 
eqs.( 2~5| , 2.15 ) a change of (f> changes the total phase of the field-operator For this reason the dynamics of <fi is 



given by the equation of motion 

• ld{H) _ 1 
9 h dN h 



An (3.2) 



where A[i is the deviation of the chemical potential from its equilibrium value. Such deviations may occur as a result 
of any fluctuations present in the system and, as discussed already, may in particular occur as a result of fluctuations 
of the value of Ao away from its average {No}. This part of the fluctuation of \i we shall denote as Ao[i- Expanding 
again to lowest order around the equilibrium No = (No) we can write 

AoM=^y(K| 2 -(AV) (3.3) 
The systematic part of the conservative part of the equation of motion of ao can now be written in the form 

( iTldi °)cons = A oM«o- (3.4) 

It is convenient to introduce the fluctuation of the free energy by AF = —TAS. The dynamics (|]J) conserves |ao| 2 
and AF. In equilibrium the right-hand side of this equation vanishes, because there A ^ = 0, and the total phase of 
the condensate (j) — (fi)t/h changes only with a rate given by the average chemical potential (/i) in equilibrium. 

The dissipative part of the equation of motion of ao near thermal equilibrium is written with the help of AF in the 
general form 

dAF(\a \2,N) 

h ( a o) dl ss = ~ r o gtf ( 3 - 5 ) 

which contains the positive phenomenological parameter To and describes the relaxation of No = \ao\ 2 to its 
equilibrium-value (No) 



According to general principles of statistical thermodynamics [h6| the relaxation process (3.5) must be accompanied 



by some form of noise. Adding a noise-term the total Langevin-equation of ao can be written in the form 

ihao = A oA1 a - iT Q dAF{ ^f' N) + f (i) exp(t# (3.6) 

aag 

Since the condensate amplitude a is a collective quantity the noise £(t) can be assumed to be Gaussian due to 
the central limit theorem. In addition we shall assume £(t) to be a white noise force. This means that the actual 
correlation time r m of the noise £ is assumed to be much smaller than the time-scale on which the dynamics of ao 
is observed, an assumption which must be checked for its validity in any concrete microscopic description. (In the 
microscopic theory we describe later it is a consistent assumption because the relaxation rate j c of |ao| 2 turns out to 
be small compared to the time-scale of motion in the trap). Thus we assume that (£(£)) = and 

(f (t)£(0)) = nk B T(2T + T')6(t) (3.7) 
(£(*)£(0)) = hk B T(T' + iT")S(t) (3.8) 



where To reappears in (3.7) because of the fluctuation-dissipation theorem. The form of the Langevin-equation 
(3.6) generalizes the work in || by taking into account a possible correlation of the phase of the condensate and 
of the Langevin force, which may exist in condensates with finite particle numbers due to gauge invariance, i.e. 
particle-number conservation. (However, it will turn out later that the coefficient L" vanishes in condensates with 
a real condensate mode i.e. without vortices, which can be understood generally as a consequence of time-reversal 
symmetry.) Gauge-invariance implies that the Langevin-equation for ao, including the fluctuating term, must be 
invariant under the transformation tf> + e. This makes it useful to write the fluctuating term as exp(icf>(t))^(t) 
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where £(t) is a complex noise source independent of </>, which, physically, describes the scattering of particles in the 
condensate with those outside .Q The coefficients V,T" describe a possible correlation of the phases of F and ato, i.e. 
the existence of a squeezing in the thermal bath of uncondensed particles, caused by the constraint of total particle- 
number conservation. We shall see that this effect actually does occur in finite condensates, i.e. the condensate mode 
imprints its (slowly) fluctuating phase on the non-condensed 'environment' due to particle number conservation in 
such a way that the lowest-lying modes are nea rly totally squeezed. 

The multiplicative nature of the noise in (3.6) raises the question in which stochastic calculus this equation should 
be interpreted: in the sense of Ito, or Stratonovich, or in some intermediate sense? This will be specified in a moment. 
Within the Gauss- and Mar koff -assumption the form of the noise-force with the same positive coefficient Lo > appear- 
ing in the dissipative part (3.5) and two further real coefficients T' , T" , is fixed by the requirement that the Langevin- 
equation must be consistent with the correct equilibrium distributio n jl 6|l P(cyq, «q) = CX P( — Ai^daol 2 , N)/k B T) 
for the condensate. Splitting into real and imaginary parts eqs . ( |3 . 7||3 . q ) b ecome 



(K(£(f))K(£(0))) = hk B T(T + T')S(t) 
(9f(£(t))9f(£(0))) = hk B TT S(t) 

(K(£(t))9ft(0))) = \hk B TT"5{t). 



Eq.( p.6[ ) may now be rewritten as 

dNo 
dt 

dt 



2 T °(N dAF 

l 



n 



k B T) 



(3.9) 
(3.10) 

(3.11) 



(3.12) 
(3.13) 



and m ust in this form be interprete d as a stochastic differential equation in the sense of Itoj]. 

Eq.(3.13) can be compared with ( |3 . 2j ) . This comparison reveals that 9?(£(t)) must describe the fluctuations of the 
chemical potential not caused by deviations of |«o| 2 from its equilibrium value but by other fluctuations in the system. 



We shall come back to this point in section VI below 



The three phenomenological coefficients Lo, L', L", are dimensionless, temperature dependent numbers, which must 
be determined from a microscopic theory. Only one of these coefficients, L , is connected with the fluctuations of 
the number of condensed atoms. If fluctuations of the chemical potential due to other processes are neglected, i.e. 
5ft£(t) = 0, the remaining two coefficients are fixed at 



r' 



r" = o 



This corresponds to the case of maximal squeezing of the noise in the phase-direction. 
Linearizing with respect to the small fluctuations 6 No <C (Nq) we find 



hSN 



2k B T 
(AN*) 



(N )To5N + 2vW^(i)) 



(3.14) 



(3.15) 



h<j>- 



d(N ) 



SN 



mm 



(3.16) 



2 A simpler ansatz (see |J) ignores the 0-dependence of the Langevin-force in (3.6). Then gauge-invariance of the Fokker-Planck 
equation which is stochastically equivalent to the Langevin equation implies F = T" = 0, i.e. the complex noise Fo(t) then has 
random phase-fluctuations which are completely uncorrelated with and equi-distributed with respect to the condensate-phase 
4>. Note however that this achieves gauge-invariance only in an ensemble sense, not for each individual stochastic physical 
realization which together form the ensemble. In contrast the form of the Langevin equation considered here does enforce 
gauge invariance for each stochastic realization. 



3 Then the Fokker-Planck equation corresponding to eqs.(3. 12 3.13) is hdP/dt = 2Y d/dN Q [N (dAF/dN Q + k B Td/dN )P} and 
has the desired equilibrium distribution ~ exp(— AF/ksT). This implies that eq.(3.6) must be interpreted in some i ntermediate 



sense, which w e ne ed not specify here further. In order to obtain its version in the sense of Ito it is best to bring ( 3. 12 ), ( 3. 13 ) 
into the form (3.6) using the Ito calculus. 
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Eq.( |3.15 ) describes the relaxation of the condensate to the equilibrium at (No) — (|ao| 2 ) on the time-scale 

h(AN^) 



(3.17) 



2T (N )k B T 

and the thermal fluctuations around it with the correlation function 

(SN (t)5N (t')) = (AA^e-l^'l/rc (318) 

The correlation time r c is an important time-scale of the problem. The noise sources Q(£(t)), 9?(£(t)) must have 
correlation times short compared to r c in order to be well described by white noise. 

On a time-scale much larger than the correlation time r c the fluctuations 5No(t) in the equation for the phase can also 
be considered as Gaussian white noise with correlation function 2T c (A.N Q 2 )5(t). Using this long-time approximation 
in the equation for the phase and taking the correlation of the effective white noise 6No(t) with Sft(£(t)) properly 
into account, <p(t) is found to satisfy the Langevin equation of a Wiener process 

d(j)(t) = y 1 D^dw 

with (dw) 2 = dt and diffusion constant 

(AJVg) d(ji) / (AJV 2 ) flfc) 1 T"\ k B T 
U * - h(N ) d(N ) { k B T d(N ) T + T Q ) + h(N ) [L ° + L) ' [6 ' W) 



l.C 



((0(t)-^(O)) 2 ) = I> |i|. (3.20) 



Eq.( 3.19 ) agrees with the result of M if we assume as in M no squeezing in the noise, i.e. T' = T" — and 
d(fj) /d(No) = k B T/ (ANq)^ 1 . Both assumptions will not be made in the present work, however (cf. also the 
corresponding discussion in the final section). 



The expectation value (ao(t)) then decays exponentially according to (ao(t)) = \J (iVo)e~ Aly l*l with the line- width 
Az/ given by the Schawlow-Townes-type formula Av = \D§- 

It is not difficult to solve eqs.( ft.l£ ,3.16) for the phase-fluctuations also on time-scales of the order of r c . The result 



for the mean square of the phase-increment in time is 

((Ht) - mf) - D,\t\ + 2^ (k B TT» + (A^)|^y) ^(e^ _ l} . (3 . 21) 
It interpolates between the diffusive long-time behavior ( |3.20|) for t > > t c and the short-time behavior for i«r c 

«*« - *<°»» 2 > = m> lT ° + T ' M + ¥ms (ms iAK) + kBTT ") '" {X22) 

The first term describes phase diffusion due to thermal fluctuations of the chemical potential on time scales much 
shorter than r c . The second term describes a non-diffusive and in principle reversible phase-collapse [||^] with the 
collapse rate 



including a contribution from the cross-correlation between both types of fluctuations. 



(3.23) 



IV. MICROSCOPIC DERIVATION OF THE LANGEVIN EQUATION 

The microscopi c der ivation of the equation of motion for the condensate amplitude «o can be carried out by using 
the Hamiltonian ( 2.29| ). As we did for the phenomenological equations in the preceding section we wish to derive here 



the microscopic equation of motion only to first order in the deviation (|ao| — y (No)) from equilibrium. 
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Hq given by eq.( 2.13| ) is the Hamiltonian, in mean field approximation, of the pure condensate. Its free equation of 
motion is 



iha = = (jm - (fj))a , 



from which 



d(f>(t) 
dt 



= -( fM) (\a (t)\ 2 )-(n))/n 



(4.1) 



(4.2) 



follows. Let us first use ( |4.2| ) in ( 2.1 8| ) to simplify H2 and then eliminate H2 by proceeding to the Heisenberg picture 
with respect to it. This changes x, x + in H3, H4 according to 



x ->■ X 



(t) = e^W J^^^e-*-"* + vta+e^) 



(4.3) 



and its adjointQ. The transformed time-dependent Hamiltonians will be denoted as H^t), H^t), but H±(t) will not 
be needed in the following. 

The equation of motion of the condensate amplitude an now takes the form,^] with the notation H(t) = Hq + 
H 3 (t) + H 4 (t), 



( dH(t) \ 
dt \d\ao\ 2 



(4.4) 



We obtain 



with 



1 dt 



h—= — Aq/x — 



( dH(t) \ 



l<*ol,X,X + 



dt 
_d\a a \ 2 

l : 

dt 



=»(?(*))- 



1 

2i' 



C/o / W(*)(e^(*)-^ (f))x(t) 



The complex noise £(t) in eqs.(|3~6-3.13) should be identified with 

£(t) =«(£'(*)) + »(&)). 



(4.5) 

(4.6) 
(4.7) 

(4.8) 
(4.9) 

(4.10) 



It is indeed independent of 4> &s required by gauge- invariance of the Langevin-equation, as can be seen from ([1.6]), 



(4.7) with ( |4.3| ).We shall see in the next section that £'(t) can be replaced by a c-number. 

To describe fluctuations around equilibrium we have replaced in the preceding expressions the quantities 
l a o| 2 ; i^Oi Xi X + by their equilibrium expressions y (No), ipo, X, X an d represent the difference in the nonequilibrium 



state by in eqs.(4.6 4.7). Omitting these differences altogether amounts to neglecting the back- 

action of the condensate on the thermal reservoir, which describes not only a modification of the fluctuating forces, 



4 For simplicity we disregard here the slow time-dependence of the f requ encies w„. 
5 For the canonically conjugate pair Nn,<f) at fixed X,X + c f- eqs.(2.( ,2.7) 
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which can indeed be neglected for fluctuations around a stable thermodynamic equilibrium, but also dissipation. 
To take the latter into account we need to calculate the averages 5{S^t{^'(t)))^, 6{5^(^(t)))ij > to lowest order in the 
interaction between the condensate and the thermal cloud of atoms. The form which these quantities must take is 
prescribed completely by the fluctuation-dissipation theorem and symmetry: 

For the reversible phase-dynamics the back-action can only lead to a shift in the average chemical potential. Such 
shifts due to the interaction H3 will be small and are neglected here. For the irreversible amplitude dynamics the 
fluctuation-dissipation theorem requires in addition the appearance of a dissipation term. If 

sjj(t-t>) = m(m(m)^ (4.11) 



is the correlation function of the fluctuating force in(4.7), the back action must modify eq.( [l.7j ) to the form 

A«o?_ mo) r di ?S M t-<)^ + 2V{No)%m (4-12) 



dt TiksT j-oo d\cto 

The derivation of this equation is given in appendix |b| 

This stochastic differential equation still differs from the phenomenological equation (3.15) in two respects 



i) The noise still has a finite correlation time r mic . We shall consider these correlation functions in more detail 
below. Taking the Markovian limit r m j c — » with 

Sjj(t - t') = hk B TT S(t - t 1 ) (4.13) 



eq.(4.12) becomes 



dt H^W^ + n^* m (4 ' 14) 



ii) The mean-field Hamiltonian i?o(|o;o| 2 ) appears in eqs.(4.12),(4.14) instead of the free energy A_F(|ao| 2 ). This 



is due to the fact that the influence of the thermal excitations on the energy are not yet taken into account. 
Doing this under isothermal or closed-system boundary conditions we should replace the energy iJo(|ao| 2 ) by 
the free energy AF(|ao| 2 ) or — TAS(\ao\ 2 ), respectively. 

This completes our derivation of the Langevin equation for the complex amplitude of the condensate. 

V. GREEN-KUBO EXPRESSIONS FOR THE TRANSPORT COEFFICIENTS 



Let us now analyse the fluctuating forces in more detail. Inserting the Bogoliubov transformation ( |4.3| ) in (^8, [1.9| ) 
the fluctuating forces take the form 

»(£(*)) = \ E (((K%T + M' v f tR ) aiala^^-^ + h.c.) 

-\-(nonre sonant terms). (5-1) 

+ (nonresonant terms). (5-2) 

Terms are called 'nonresonant' if the frequencies of the quasi-particles cannot add up to zero. Such terms have not 
been written out explicitly because later-on we shall restrict ourselves to the resonance or rotating wave approximation 
in which they don't contribute. The relevant matrix elements M^\M^ are 

M { K % = 2U J d 3 x^«^ + + {u <-> /*) 

(5.3) 
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and very similarly 



(5.4) 



The matrix-elements M'W, M'( 2 ) coincide with M^ 2 ' if the dependence of tpo on (No) is negligible or vanishes, 

as e.g. in homogeneous systems. 

Mj^lfj, and similarly M'J-^ describes a scattering process in which one atom is scattered out of the condensate by 
the absorption of the two quasiparticles u, \l from - and the emission of the new quasiparticle k into - the thermal 

(2) '(2) 

bath. Likewise M^' K and similarly M^' K describes a scattering process where an incoming thermal quasiparticle 
k is absorbed, again an atom is kicked out from the condensate, and two quasiparticles v, fi arc emitted into the 
thermal bath. The scattering amplitudes for both processes are linearly superposed due to the phase-coherence of 
the condensate which exists on the time-scale of the relaxation process induced by the scattering process even if it is 
destroyed on a much longer time scale. 

We can now calculate the correlation functions of the fluctuating forces. Their averages over the bath of quasi- 
particles vanish, (3fi(£(i))) = = (3(£(t))). Their second-order correlation functions are obtained as 

mi'itmi'it')))^ = 1 e \{(m' k % r +m$ k )\ 2 

U K (n v + l)(n M + iy("~-^-^)(*-0 

+n v n„{n K + i) e -^-^-^){t-t') j (5 5) 
|n K (n, + l)(n„ + l) e ^-^-<^)(*-*') 

+n„n M (n K + i) e -*(o»-o»-^)(*-*') | (5.6) 

<3i(£'(*))3(£(0)>* = ^ E { (i M ^%T ~ M Hn) (<% + «%)*) 

(n„ + ljfa, + i)^ e ^-^-^)(t-0 

These correlation functions can be replaced by delta-functions provided that the frequency-sums contain a flat quasi- 
continuum of nearly resonant terms in a neighborhood of the resonance u) K — Cj v — u)^ — which is broad compared to 
the damping rates we calculate here. This assumption will be satisfied in s uff icien tl y lar ge condensates. The strengths 



(5.7) 



of the delta- functions can then be extracted from the expressions ( |5.5| ), (p.q), (5.7) by taking the time- averages 

L°L d(t - tf)(i(t)i(t'))^ and d(t - n(i + (tW))<f>- 

•Ji(£'(i)), 3(£(i)) are here given as expressions involving operators. Provided the Markovian approximation is satisfied 
the average of their commutators over the quasi-particle bath are again given by delta functions in time. Explicitly 
we obtain for the coefficients of the delta-functions 
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d«f(i)),sR(r(o))]> = o= / dt(m(t))Mm)})4> 

OC J —OQ 

{(nj, + !)(% + l)n K - (n K + l)n„n, M }5(5;« - u>„ - u^) = 



(5. 



It can easily be verified that the bracket {...} in the last line of (5.8) vanishes if it is multiplied by the delta- function 
expressing energy-conservation. As a result the fluctuating force £ in the Markoffian limit can indeed be treated as 
a c-number and will henceforth again be denoted by £. This also serves as a nice consistency-check that it is indeed 
possible to treat the condensate classically, even after taking its interaction with the quasi-particles into account. 
Let us now proceed to derive formulas for the three transport parameters To,! 1 ' and V" . From 



2r + r' = 
r' + it" = 



1 



fiksT 
1 

hksT 



— oo 
+oo 



dtmm)* 



implied by eqs.(3.7,3.8) we obtain the representations 

1 



To 

r + r' 

r" = 



1 

TiksT 
2 

TiksT 



OC 

oo 



OC 
OC 



d(t-t')(mt)mat'))u 

d(t-t)(St(?(t)M?(jf)))4> 

d(t-t')m'(mim)u- 



(5.9) 

(5.10) 
(5.11) 
(5.12) 



which have the form of Green-Kubo relations for the transport coefficients. Using the explicit fo rm ( pM|]5.2| ) of the 
fluctuating forces the thermal averages can be taken and the time-integrals in eqs. can be carried out 

which leads to the formulas 



r ° = 2^T E \( M i%T - MW K \ 2 n„n^n K + 1)S(Q K - £7, - w„) 



(5.13) 
(5.14) 



r" 



2hk B T 



E ((Mi%T-M% K ) (M' K % + (M^J 
^n M (n K + l)S(u K -tip- ui v ) 



(5.15) 



These expressions constitute our general results for the three transport parameters. They have to be evaluated 
separately for each individual trap geometry. 



VI. RELATION TO THE FLUCTUATION AND DISSIPATION OF THE EXCITATIONS 



As was pointed out after eq.( 3.13| ) by phenomenological arguments, the noise term 5ft(£(i)) is not connected with 
the fluctuations of the number of particles in the condensate, but must be due to other fluctuations, which are then 
necessarily thermal fluctuations of the amplitudes of the excited states. In our microscopic results this can be seen 
from the fact that the fluctuating force !ft (£'(£)) according to eq.(L8) contains precisely the same operator which also 



appears in H 3 (t) and couples the atoms in the thermal cloud to the condensate. 

In the special case where the difference betw een t he coupling matrix-elements M'^ 1 ' 2 ^ and M^ 1,2 - 1 is negligible (which 
is exactly satisfied in box-like traps, cf.section VII) the intensity T + V of the noise source can be expressed entirely 
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as a property of the excitations, as we shall now demonstrate.^ For their amplitudes a u (t), 6t^(i) a quantum-Langevin 
equation can be derived microscopically along the same lines employed here for the condensate amplitude. We have 
done this elsewhere p^j, see also || with the result, in the Markoffian limit, 

= -iw v a v {t) - j v a v {t) + |„(t) (6.1) 
with Gaussian fluctuating force-operators with vanishing average and 

(6.2) 

([6(^(*')]>=27^(*-*')<W 

where the damping rates 7„ are given by 

> = ^ E {\( M i%r + J a (n M - n K )<5( WK - - w „) 

+ (^,)1 2 (n« + + c M - a;,)}. (6.3) 

The first term describes Landau-damping of the mode v by scattering a quasi-particle from mode fj, to mode k and 
is equivalent to a result derived in Jl9| by the golden rule. The second term in eq. fl6.3| ) describes Beliaev damping, 
where the mode v decays into two modes ft, [i. It survives even for T — > where n K — > for all mod es. 

Let us now establish the connection between To + T' and the damping rates y„ as given by (3.3). We shall show 
that the simple sum-rule 

ro + r-^A^^^ + ih, (6.4) 

holds. To see this we need to consider 

£fH,(n„ + 1)7, = ^> £ |(M#„)* + ^.J^ - w - - ^) 

■{(%, - n K )n„{n v + 1) + ^(n M + n„ + l)n«(ra« + 1)| (6.5) 

The second term in the curly bracket arises from the second term in (3.3) by first exchanging the notations for the 
summation indices v and n and then symmetrizing in v and /i, because the matrix-elements are already symmetric 
in these indices. The remainder of the proof then consists simply in noting that for u> K = u) v + lo^ the identities 

- n K )n v (n v + 1)= n li n v (n K + 1) 

(6.6) 

(ra p + n y + l)n K (n K + 1)= n li n v (n K + 1) 

hold. Using this in eq.(3.5) and then comparing with (|5.14 ) establishes the sum-rule. We can also note that the 
processes due to Landau scattering contribute to the sum-rule with precisely twice the strength of those due to 
Beliaev scattering. 

Thus we see that, in general, the noise- amplitudes proportional to the combination of matrix-elements {M^l^)* — 

Muf) }K drive the ratimfeer-fluctuations in the condensate, while those proportional to (Af*,^)* + Mvp, K are due to 
fluctuations of the occupation numbers in the excited states, couple in the Hamiltonian to the particle number in the 
condensate and therefore drive the p/iase-fluctuations in the condensate. 



In the general case the coupling of the condensate to the non-condensate modes differs from the coupling between the 
non-condensate modes and the relation between To + V and the -y v is less direct. 
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VII. EVALUATION OF THE TRANSPORT PARAMETERS FOR A BOX-LIKE TRAP 



For simplicity we consider now a trap consisting of a cube of length L with cyclic boundary conditions. In the 
following equilibrium values of all parameters are implied, but we shall omit in this section the tilde and write /j, for 
(/i) to simplify our notation. The normalized u and v coefficients are in this case 



u v = 



E v +pl/2m 1 i^.s/n 



V y/2E v pt/m W 



E v - pl/2m 1 
y/2Ejpi/m W 



Jpv£/h 



with 



and p v — h^n u with integer vector n v . 



(7.1) 



(7.2) 



(7.3) 



A. Transport coefficients 

The squares of the relevant matrix elements for E K — E v + E^ become 



Ml 2 .) 



(m« Y- 
(mW Y + mW 



V J E V E^E K 



UoV G(E V7 E M - f i) 
V J E V E^E K 



with 



G(x, y, a) = \J a 2 + (x + y) 2 ^3V a 2 + x 2 \/a 2 + y 2 — xy + 2a{\] a 2 + x 2 + \J a 2 + y 2 + a)j 
+ (x + y) (x{\J a 2 + y 2 + a) + y{y/ a 2 + x 2 + a)j 



+2a(y / a 2 + x 2 + a)(y/ 'a 2 + y 2 + a) + a(x 2 + y 2 + a 2 ) 
The transport coefficients are then expressed as the simple result 

r" = o 



and 



2 (a\ 2 ^ 8{ e v + e M - e v+li ) t 



r -;(z) E 

r + r'= - (-)" ^2 6i: " :,/+,,) 



k T 

G{e v , £ M , a)F(e v , e M , - — ) 



k T 

G\£ v , £u, — Cx)F[£ I/ , £„, — ) 



with 



PC k BT , 



e /3hoj (£„+£,») 



TiujQ k B T [e< 3huo ^"+ e ^ - l) (e^o £ " - 1) (eP nMoe » - 1) 



(7.4) 
(7.5) 



(7.6) 



(7.7) 

(7.8) 
(7.9) 

(7.10) 



Here we scaled the scattering length a = mUo/4-Kh 2 with L and the energies E^^E^ and (j,, ksT with the energy 
huj = {2ttK) 2 /2mL 2 , defining 
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e v = yj (n 2 + a) 2 - a 2 



(7.11) 



e„+„ = y / ((^ + ^) 2 + «) 2 -« 2 ( 7 - 12 ) 

with a = fi/hujQ. 

The double sums over n^n^ start with n- values with |n| = 1. They are approximated by integrals according to 

OO OO 

V 5 ( e „ + s M -S„ +M ) (...)= tt 2 / / £,£„(£„ + £„)(■■■)&„&„ (713) 

yrWyit V /(4 + « 2 )( £ 2 I + « 2 )((^ + ^) 2 + « 2 ) 

Here (...) is any smooth function of e v ,e^. In all experiments so far a > 1 is satisfied, i.e. we can replace 



yl + 2a — > y2a. This leaves us with the integral expressions 



OO OO 



a\ 2 f f G(e v ,s li ,a)F(e v ,s li ,^)de 1/ de M 



p Q=27r - , , - - ™' - = (7.14) 

' * i- J(el + a 2 )(el + a 2 )((e„+eJ 2 +a 2 ) 

2a v 2a V 



OO OO 



- i- ,/(£2 + a 2 ){el + a 2 )((e v + £ M ) 2 + a 2 ) 

a V2a V 



The expression for To and the asymptotic behavior for 6^,5^ — > 0: G(e v , e^, \a\)F(e v , e M , hsT/TUPo) —> 
18(kBT/hu}o) 2 l/[e u e fl (e lJ + £ M )] make it amply clear that the states with the smallest energies £„ <C a make a 
large contribution to T (but no£ to r + V). To calculate this contribution it is permitted [] to use Phcooe^ <C 1, 
(3TiojqE^ <C 1 under the integral and to approximate 



knT , / fcflT N 2 



F( £ „^,^-)= -f- 7 — — (7.16) 

ftw \ "-Wo / (e„ + £ M )£„£ M 



in addition to approximating \J e 2 + a 2 ^je 2 ^ + a 2 yj (e v + £ M ) 2 + a 2 = a 3 , and neglecting terms of order e 2 /a 2 . 
This contribution to To, which we shall denote as Too, then reduces to 



r o = 36^) 2 (^n / de v [ de^ J- (7.17) 

\LJ \huj a J J J + Eft) 

The double- integral can be evaluated as \j2ja In 2 which yields the final result 

r 00 . & (*) » (£)' . , 5 , g) 2 )"V «, 18) 

The second form of this expression is obtained by eliminating V = L 3 in favor of the critical temperature of the 
equivalent ideal Bose-gas at the same density via 

V = (N/({3/2))(2irh 2 /k B T c m) 3 / 2 . (7.19) 

There is yet another particularly important contribution to To due to the infrared-singularity of the integrand, we 
shall denote it as Toi, where only one of the two excitation frequencies, say E v , is small compared to /i, while the 
other is larger, of order fj, or even ksT. With respect to e„ the low-energy asymptotics may then still be used. We 
then obtain 



7 We actually need the additional condition \J JiEuJo <C ksT 
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r i 



47r / a \ 2 
a 



qa 



de v I de 



G(0,e M ,a) - G(0,0,a) e ^° e * 



(7.20) 



.'2« 



/2a 



where we have set an upper cut-off for the small energy at a fraction q of the chemical potential. Most of the e^-integral 
comes from a range around a and we may therefore replace the thermal function by its asymptotics for phuoe^ — ► 0, 
which is (k B T /TvjjQe^) 2 . The integrals can then be performed with the result 



01 (l) i y 2fiu!o} hujo/j, 



1.22. 



T\ 2 N (k B T c a 2 m 



T, 



(No) 



1/2 



In 



2Hujq 



(7.21) 



We conclude that this contribution to To is smaller than the leading term Too by the order of magnitude 
(hu} /ti) 1/2 \n(q 2 n/2hLJo). 

Let us now turn to the expressions for T — T 00 — r i and r + V. They can be simplified for a ^S> 1 by rescaling 
£„, by a and taking the limit \]2j a — ► for the lower boundaries of the rescaled integrals. We find in this way 



ro-r 00 -r 01 = 2 7 r(|) 2 (£-) J dx J dy 





G(x,y,l)F(x,y,^f) 
^ (x 2 + 1){y 2 + mx + y) 2 + 1) 

k B T\2/ 18 



+8- 



fi / \xy(x+y) 
-\x 2 + l)(y 2 + y/y 2 + 1 - 1) + y{y 2 + l)(x 2 + Vx 2 + 1 - 1) 



x 2 y 2 {x 2 + l){y 2 + 1) 



^ L J ( Twoq 



dx I dy 



G(x,y,-l)F(x,y 



^(x 2 + l)(y 2 + i)((x + y) 2 + l) 



(7.22) 



(7.23) 



where we used the scaling property 



„, k B T 1 e v k B T x 

Uluo a a a nuj^a 



(7.24) 



Eqs( 7.22 , 7.23| ) are the complete result for the temperature dependent transport parameters of the condensate for 
box-like traps. In general the integrals have to be done numerically. We shall here consider some asymptotic results 
only. 

First we consider these expressions asymptotically for k B T / /i 3> 1. Then the integrals receive important contribu- 
tions from x, y of the order of 1, i.e. from quasi-particle-energies of the order of the chemical potential, and also from 
values of x,y large compared to 1, i.e. quasi-particle energies of order k B T. The contributions and r'( > ) from 
large energies can be determined in leading power in (k B T j \x) by approximating 



F(x,y 



k B T x 



-(x+y) 



k B T 



(7.25) 



and rescaling x and y by k B T//i. In the integrals^] for To and To + V we can then let fi/k B T — > without any 
problem, using the property G(x, y, 0) = 4xy(x + y). whereupon they are easily evaluated with the asymptotic results 



r (» + r '(»^ r (» 



*(z) 



a\2 k u T 



Tlujq 



1.27. 



T k B T c a 2 m 



(7.26) 



We see that the result for vanishes to this order. 

For the contributions Tq ,T'^ from quasi-particles with energies around \i we can approximate F(x,y,k B T/iJi) 
according to eq.(7.16) and find 



In the high-energy regime the subtractions of the infrared-divergent terms in the integrand of (7.22) are of no importance 
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r <^ B w ( .f t^ii _ . 0243 .. s <„, (ry^ {T21) 



*LJ Tlluoh ' " u \T C J (No) V ti 
(u) >(,a ,(,,W a\ 2 (k B T) 2 ,,„, / T\ 2 N ( k B T c a 2 m 

VL/ hcjofi \T C J (No) V ?i 2 

with the numbers Bq , defined by the integrals 



1/2 



(7.28) 



o o 



^/(z 2 + l)(y 2 + 1)(0 + y) 2 + l>y(x + y) 
18 

a;?/(x + y) 

Q a;(a; 2 + l)(y 2 + V?/ 2 + 1 - 1) + Z/(z/ 2 + 1)(^ 2 + V^TI - 1) 



(7.29) 



x 2 y 2 (x 2 + l)(y 2 + 1) 

oo oo 

B'M =2n [dx [ dy - G(x,y,-1) 

7 7 A/(a: 2 + l)(w 2 + l)((x + w) 2 + l)a;w(x + y) 



o o 



We can conclude that the contribution from quasi-particles at energies of order fi is larger (for To by an order 
of magnitude k B T j \i) than the contribution from energies of order k B T, but Iq is, in large condensates, still 
subdominant to Too by the order of magnitude \Jfuooj \i. 

1/2 

Now let us consider also the low-temperature limit, namely k B T '/ fx <C 1 or, equivalently, T/T c <C (k B T c a 2 m/T% ) 
In this region it is not neccessary to distinguish N and (No). The inte grals now receive their contributions for x, y 
both small compared to 1, but we can still use the approximation ( 7.25| ). For small x,y we can expand 



G(x, yi l)~l8 + 3(x 2 +y 2 +xy) G(x,y,-l)~^x 2 y 2 (x + yf (7.31) 

To obtain the leading term it is enough to keep only the smallest powers of x,y in the integrands. The integrals are 
easily evaluated with the asymptotic low-temperature results 

„ „ / a\ 2 k B T T knTra m 

r.-r M = 3to(j) ^ = 5.72..-^- (7.32) 



As long as the temperature is high enough to satisfy k B T 3> y/JIhujo the part Too still dominates the value of To 

B. Particle- number fluctuations 

We follow the procedure of Giorgini et al. [Q and deduce the particle-number fluctuations in the condensate from 



the number fluctuation in the thermal cloud. This leads to eq.( 2.28] ), which we evaluate using the expressions for 
u v , v v and E v . 
We obtain 

(AN 2 ) = £ 2M ^ + 1)( ^ 2 +V)+m2 (7.34) 
Approximated by an integral this becomes 



de i \Je 2 + a 2 - a ( 2 s 2 + 2a 2 
^i+^TV e 2 + a 2 [ a + 2(sin/i^2£^- 



(AN 2 )=k -i - ^ , ^ a< + — -gg— (7.35) 
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The dominant contribution comes from the lower boundary of the integration |L4j which contributes, for a»l 

(A iV 2 ) ~ 2n (^j = A' ( Hlhf-Xv^ (7.36) 

with 



h 2 



A' = — ^ = 0.0161.. (7.37) 
Ztt 

More precisely the dominant contribution to (AJVq) is given by the discrete sum [E3 

<A7V 2 )=2 Ai 2 (fc B T) 2 ^i T (7.38) 

which gives the same expression as ( 7.36| ) but with the prefactoiQ 



A 2 



J2 -4=0.021.. (7.39) 



(2tt) 4 ^ nf 



If we eliminate the volume in favour of the critical temperature of the ideal Bose-gas of the same density via eq.( 7.1S ) 
we get 



C(3/2) 4 /3 y T( 



At temperature T = a similar evaluation of ( 2.28 ) gives 



(A7V 2 )| T=0 - 2^{aNf' 2 V~ 1 / 2 . (7.41) 



C. Particle- number relaxation rate 



We are now in a position to evaluate the rate t c 1 from eq.( |3.17| ) using the results for {ANq} (numbers are calculated 
with the prefactor A') and To ~ Too- We obtain 



This result applies also in the low-temperature region, because it makes use only of the results for (ANq) and To 
which also hold in that region. 

To get an idea of order of magnitudes we compare this and the following results with the damping rate 70 of the 
lowest lying modes, which is given by |l9| , |2"l"| ] 

70 = ^ ( 1 ) **I = 4.06..^-i/3 ( ksTWrrA ^ M 

We see that the relaxation rate j c is of the order 

»ri f kBT c a?m\ . , 

7c „ iVi/3 I JLj. \ 7Q (7.44) 



The proportionality factor is of the order of •v/ h/Tujo and is large in large and strongly interacting condensates. I.e. 
then the relaxation of the condensate to its equilibrium is faster than the relaxation of the low-lying collective modes, 
but slower than the frequency of the lowest lying modes, which is y^2uJo/j,/h. 



9 The formula ( fr.39 | ) di ffers from the one given in Jl4| by a factor 2 4 whereas the numerical result differs by yet another 
factor; the formula ( |7.4l| ) differs from the one in [jl4| by a factor 2. 
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D. Phase collapse rate 



The phase collapse rate is given by eq.(3.23) and requires only the result for (A Nq) and T" = 0. At zero temperature 
it reduces to 

IcollapseW^ = ^^)V< AJV 0>lT=0 (7-45) 



i 23.6. v 3 /4^a k B T c f k B T c a 2 m\ 7/S 1 
7c ;jap Se |T=o = — ==-(ara ) 7 — = 2.50..—— -= (7.46) 



from which we get 

y/V m h \ H 2 ) y/N 

For finite temperature we obtain 

^bT a , l nfk B T c a 2 m\ 1 / 2 



OQ7R k B l 1/3 /fe B l c a-m\±/^ 
Icollapse = 0.876..— — N > * y J (7.47) 



By comparison with ( 7.43 ) we see that ^collapse is of the order of 70 and is therefore in large condensates smaller than 
7c- 

In summary the phase-collapse is not effective in large condensates because it occurs with a rate ^collapse < 7c and 
is at the same time restricted to a time interval At < — since for larger times phase-diffusion takes over. 



E. Phase-diffusion rate 



The phase-diffusion coefficient is a somewhat complicated quantity because it receives contributions from several 
processes, which are physically distinct. We consider the different contributions separately and also distinguish the two 
temperature regimes of high temperature k B T > fj,, for which we give the result first, and k B T < fx (low temperature). 

1. Low frequency condensate number fluctuations 



From eq.(3.1£) we infer with Tq = T. 



no 



2 



which is evaluated as 



D-= 0.0853..- J^p-j ^^TT^iT- < 7 - 49 > 



The same result holds in the low temperature regime k B T < /i. In comparison with 70 (7.43) it is of the order 



70 ~W— 7o (7.50) 



fc ) f /1 

and is much smaller in large and strongly interacting condensates. Still this contribution to the phase-diffusion rate 
is always the dominant one at low temperatures and may dominate even at higher temperatures (see below). 

2. Condensate number fluctuations due to quasi-particles around energies /j, 
Splitting To = Too + (To — Too) and expanding to first order 

1 1 r - r 00 



To Too T 



(7.51) 

00 
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by (B 



we estimate as contribution from the higher frequency condensate number fluctuations described by Tq — Too as 
given 

D * r e <JV ) » ^ 7 - 52 ) 



which is in absolute value smaller than the contribution from low-energy excitations by the order of magnitude 

factor yj Tiluq / [i. This contribution is therefore negligible in very large condensates. In not so large condensates the 
complete integral in the result for Tq needs to be evaluated. 
In the low-temperature regime fc^T < fj, we get instead 

<)_o.306..i(M>!^y /2 M (7 53) 

This is much smaller than by the order of magnitude factor (^i/fe^Te) 1 / 2 /TV 1 / 3 . 



3. Fluctuations in the thermal cloud at energies of order fi 



By eq.(3.19) this contribution is given by 



which is evaluated to 



M ^/T^ 3 N fk B T c a 2 m\ l2 k B T c 

* = °- o243 -^ ) irj w r^H ^r- (7 - 55) 



This contribution differs from by the order of magnitude factor {T/T c ) y/JjMUo/kBT c and is therefore much 

smaller. 

For temperatures k B T </iwe find instead 

which is again negligibly small compared to D^"\ 

In summary, the phase-diffusion is caused dominantly by the lo w-frequency particle number fluctuations in the 
condensate and the phase-diffusion constant is given by eq.(7.49). It is proportional to temperature, and scales 
proportional to N~ 2 / 3 for fixed T c or proportional to iV -1 / 2 for fixed volume of the trap. 



VIII. EVALUATION OF THE TRANSPORT PARAMETERS FOR AN ISOTROPIC HARMONIC TRAP 

In this section we consider the more realistic case of condensates in a parabolic trapping potential itiluqX 2 /2, which 
we assume to be isotropic for simplicity. In order to analyse the noise 5(£(f) driving the fluctuations of |oo | 2 we must 
consider in detail the relevant linear combination of matrix elements 

(Mi%)* - = 2U J d 3 x^ {(u K - «„)(«*(;* + - n^v* + v R v;v* v } (8.1) 

In the following we shall make use of the local density and Thomas-Fermi approximation, restricting ourselves to 
large condensates. For the high-lying states we can then use the local energies in Thomas-Fermi approximation 



E(p, x ) = \l^ + \Uon (x)\) 2 ~ U 2 nl(x)Q(p - V(x)) (8.2) 

with the condensate density 
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n (x) = (iV )|^o(^)| 2 = ({l4/Uo)(l - {x/r TF f) (8.3) 

and the Thomas- Fermi radius 

rr^^^^i 1 -^^)^. (8.4) 
The high-lying quasi-particle modes can be represented similarly to the spatially homogeneous case as 



c/n 



(8.5) 



_E ! L -p^/2m plpK . x/n 
y/2E K p*/m 



The low-lying collective modes can be represented as 



u Q n Q (x) i / m v , 



(8.6) 



. U n Q (x) 1 / M) v . 



with 



d A x\ X v{x)\ z = 1 (8.7) 



p3p 
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form 


m 



Xu{x) = -^p£ n " ] {x/r TF )(x/r TF ) l Y lm {9,ip)Q{l - x/r TF ) (8.8) 

r TF 

The polynomials Pj- 2 ™' 1 (x) of degree 2n are the normalized solutions of the radial part of the Bogoliubov- Fetter 
equations in the Thomas-Fermi and long- wavelength limit J23|,0 given by |24 



p? n \*) = V4n + J £ + 3 -- 2e ~ 1 ^ [* 2n+2e+1 ft - * 2 r] (8-9) 

with the normalization 

f dxx 2e+2 [p( 2n \x)} 2 = 1 (8.10) 
Jo 

~— 1/2 

In the phonon part of the excitation spectrum we have u\ ~ — wa ~ w A . Furthermore, in that low- energy 
region the statistical factor in eqs.( |5.13' ), (5.14), Q5.15 ) is well approximated by n^fi^n^ (ksT) 3 jTi~ H) k lo v lo M . Just 
as in the case of box-like traps the frequency factors in the denominator, together with similar f urthe r factors in the 
denominator coming from the matrix elements, make the phonon contribution to the sums in J5.13j ) the dominant 
one, at least in large condensates, and we shall therefore concentrate on this contribution in the following. This 
frequency range has a natural upper cut-off at (/i) /%, where the collective phonons go over smoothly into particle-like 
excitations. 

For E K , E U ,E^ -C (/i) the matrix elements (M^)* , M^f} tK are given by the integral 
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(Mi%r * -m^ii ( 8 - n ) 

V 8vr r \ F ^/2E v E ll (E v + E^) 
where J denotes the integral 

JK,„„,^;4,4,^) = f dxx 2 {l - x 2 fx l ^ l ^P^\x)P^\x)P^\x) 

Jo " 

and the Clebsch-Gordon coefficients C'(l K |£ Jlt m /J , l u m v ") are given by the angle-integral 



if \£fj_ — l v \ <i K <t^+lui otherwise they vanish. Later-on we shall have to calculate e.g. J2 m „ m „ K-^'V'm)* — M^ k \ 2 
where we can make use of the sum- rule, for \£^ — 1„\ < l K < + t v 

W^m^m^f = 1 

so that the Clebsch-Gordon coefficients need actually not be used explicitly. 

In order to have well-defined expressions for the rate-coefficients we again need to smoothen the delta-function 
expressing energy-conservation, which is done physically by experimental imperfections or limitations in resolution. 
Here this can be done by replacing the discrete sum over the 'quantum-number' £ K by an integral 



dE K jdl K 

^ ^^ ^^-W (8 - 13) 



where we used the expression for the excitation-energies ]23| 

E K = hu> e K 



e« = y/2n 2 K + 2nJ K + 3n K +e K . (8.14) 

We introduced the dimensionless eigenvalues e KjIAA1 which will appear in the ensuing expressions from now on. The 
integration over E K with the delta-function then picks out the energy-value E K — E v + E^ so that t K becomes a 

function Pk of the other 'quantum-numbers' 



{e v + e^) 2 - 2n\ - 3n K 



2n K + l 
The inequalities 

then imply that n K must lie in the interval 

n K - < n K < n K+ 

with 

n K ± = ^(yj \l v ± l»\ 2 + \t v ± t»\ + 9/4 + 2(e„ + e M ) 2 - % ± ^| - 3/2) 
Using all this we obtain from eq.( |5.13 ) 



„ „ / a knT \ 2 T 2 kBT r a 2 m , N x 1/3 , 

r ° = M^) = Boo ?? J SMc(^ (8 - 15) 
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where Tiloq is now eliminated in favour of fcgT c via the relation Hluq = fcsT c (C(3)/A^) 1 / 3 , and where the temperature- 
and particle-number-independent positive real number Bqo is defined by the multiple sums 

B ™-—l^l^l^l^ 2. elel(e„ + e,)(2n K + l) (8 ' 16) 

The result for To agrees, except for the numerical prefactor, with the result of Q which was evaluated there using 
the local-density approximation and imposing a lower cut-off for the excita tion- frequencies at the geometrical mean 
trap-frequency Q. It can also be compared with the corresponding result ( (7.1S| ) for the box-like trap, which shows 
the same dependence on temperature and particle- number (if we stipulate (No) ~ N), but the comparison of the 
prefactor is problematic because the condensate in the parabolic trap has two length-scales do and ttf, whereas in 
the box- like trap only the length-scale L is relevant. 

The property flS.ll ) of the matrix-elements implies that the low-lying excitations don't contribute to Tq +V. The 



reality of the matrix-elements furthermore implies that T" vanishes. These remarkably simple results mean that the 



noise-source £(t) introduced in eq.(3.6) is purely imaginary, corresponding to total squeezing in the direction of the 
phase 4>. In other words, the coupling of the condensate to the collective excitations introduces a direct Langevin 
noise-source only for the number-fluctuations SNq, not the phase- variable 4>f^\ 

The fact that ro + V = for the contribution from the low-lying states implies that the contributions from the 
higher lying states must also be considered in order to evaluate the small (compared to Tq) but finite value of this 
quantity. For this purpose we need to consider the matrix-element (M'^J^)* + M'J 2 ^ K .. It differs from the matrix- 
elements we have considered so far by the replacement ipo(x) — > ipo(x) + 2(No)dipo(x)/d(No) in the matrix- element. 
In Thomas-Fermi approximation this is tantamount to the replacement 

Mx) - (2/5)^0 (x)/(l - x 2 /r 2 TF ). (8.17) 

Physically this implies a reduced coupling of the thermal fluctuations with the center of the condensate and a 
strongly enhanced coupling at its boundary, as one would expect for fluctuations located in the thermal cloud. A 
mathematical consequence is the fact that the integrals defining these matrix-elements diverge at the boundary in 
the Thomas-Fermi approximation, meaning that we encounter here the limitations of that approximation. Instead of 
a full-fledged extension of the theory beyond the Thomas- Fermi approximation it will be sufficient for our purposes 
here to cure its deficiencies by substituting as a cut-off the finite thickness of the boundary-layer given by 122] 



/ huJo_ 



d = —ttf 



2/3 



The matrix-element itself is then evaluated in the local-density approximation |22] , where we can make use to good 
purpose of the analysis already performed in the predeeding section. The finte volume V = L 3, (and the associated 
Hluq = (2irh) 2 /2mL 2 which is not to be confused with the trap frequency called ojq in the present section) is then 
an arbitrary local subvolume of the condensate, introduced merely as a technical device like a quantization-volume. 
It must be sufficiently small so that the condensate within it can be treated as homogeneous, and sufficiently large 
that we can replace sums over local momenta by integrals. At the end we have to check for consistency whether 
the result is indeed independent of the choice of this volume. The result obtained in this way is the local average 



of the result (7.28) for the homogneous case, which now becomes space-dependent, because we have to substitute a 
space-dependent chemical potential /i — > — x 2 /r' TF ). This local result can be written as 

rM ^r'M B ' ( " ] (k B T) 2 a 2 m 
To(x) + T (x) - 



2tt 2 a 2 < M )(i-z 2 /4F) 



and is indeed independent of the choice of V. The local average has to be performed with the weight (ipo(x) + 
2(N )d^ (x)/d(N )) 2 determined from fcvj\) . Doing the average and regulating the divergency of the integral at the 
boundary of the condensate by the physical cut-off we obtain 



10 The latter is of course affected by the noise-source indirectly, because the fluctuations of 5iVo driven by the latter cause 
fluctuations in the chemical potential. 



23 



r n + r' 



2 l/3 
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hujQ J 



(No) 



-2/15 



28/15 



0.024.. B l{tl) 



N 

W) 



2/15 



N 2/9 



k B T c a 2 m\ 



14/15 



(8.18) 



In order to extract results for the relaxation-rate of the condensate number and the phase-diffusion rat e it i s 
necessary to know also the mean square of the number-fluctuations (ANq). This can be evaluated from eq.(2.28), 
using the fact that these fluctuations are also dominated by the low- lying modes 14 1. The result of this calculation 
to leading order in (haj Q /k B T) is 



A 



T 

(C(3)) 8 /i5 \T C 
with the number A given by the multiple sums 

2£ + l 



(No) 
N 



4/5 



AT4/3 



e B l c a m 



2/5 



n n' 



(e(n,£)e(n',e)y 



dx(l-x 2 )x 2 ^Pf n \x)Pf n '\x) 



(8.19) 
(8.20) 

(8.21) 



In order to find the scaling of (ANq) in the thermodynamic limit N — > oo, ujo — > 0, k B T c = huJo(N/ ((3)) 1 / 3 fixed it is 
necessary to use the form of the preceding results in which hu>o is eliminated in favor of k B T c and to use (No) ~ N.. 
Then the scali ng (A N 2 ) ~ N A ^ derived in is recovered. The particle-number relaxation rate now follows from 
eqs.( 3.17 ) and ( 8.15 ) as 



7c 



2B 
A 



■(No 



.1/6 (JL 

d Q 



6/5 k B T 2(C(3)) 1 /5 jBoo T /<iV )\ 1/5 /fc s T c a 2 m\ 3/5 k B T c 
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n 
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It is the largest of the various rates we calculate here but is still small compared to ujq, the inverse time-scale of motion 
in the trap, by the order of magnitude N~ 2 / 3 (N a/do) 6 / 5 . 

The phase-collapse rate is obtained from ( |3.23| ) . At T ^ ( more precisely above a cross-over temperature of order 
hu>o) we find 



"/collapse 



15 V5 A i/2 5 ( ± y 6 kBT = 15 2 / 5 (C(3)) 2 / 15 AV 2 T (_N_y /B N -i /3 ( k B T c a 2 m\ 2/5 k B T c 



d 



T c \(N ) 



.23) 



Apart from the numerical prefactor this is the same asymptotic expression as obtained for the damping-rate 70 of 
the low-lying collective modes (see e.g. It is smaller than j c by the order of magnitude ((No)a/do)~ 2 ^ 5 , i.e. the 

phase-collapse remains inefficient before phase-diffusion takes over. 

The phase-diffusion constant D^f 1 due to the exchange of particles between the condensate and low-lying excitations 
is gotten by inserting the results for (AN 2 ) and Too in eq.( 3.19 ): 



25B 
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It is smaller than "/collapse, again by the order of magnitude of ((No)a/do) 2/5 . 

Finally, the contribut ion o f the fluctuations in the thermal cloud to the phase-diffusion is also obtained from ( 3.19| ) 
by inserting the result flS.18|) for T + T': 
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(8.25) 



It differs from the previous rates, which were all proportional to temperature, by the stronger temperature- 
dependence ~ T 3 . However, this contribution to remains smaller than by the order of magnitude 
W- 1 /9(fc B T c a 2 m/?i 2 ) 11 /i5(r/r c ) 2 . 
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IX. DISCUSSION AND CONCLUSION 



In this paper we have put forward a detailed theory of fluctuations and relaxation processes of the condensate in 
thermal equilibrium with the cloud of its excitations. For a given number of particles No in the condensate, we have 
defined the condensate mode as the corresponding normalized solution of the Gross-Pitaevskii equation, defining at 
the same stroke the -/V -dependent part of the chemical potential. The equilibrium value of (No) is distinguished as 
the value of No for which the number of particles in the thermal cloud in equilibrium with the condensate plus No 
equals N. We have calculated the fluctuations of No around its equilibrium value and also the fluctuations of the 
phase of the complex amplitude ao of the condensate with \ao\ 2 = No- In a general phenomcnological framework 
presented in the first part of this paper we were able to separate the fluctuations of the complex condensate amplitude 
into several contributions, which have different physical origin: 

- The fluctuation of the atom-number in the condensate, which are driven by the exchange of atoms between the 
condensate and the thermal cloud. 

- The fluctuation of the chemical potential with two different contributions, namely the fluctuations of \i due to 
number-fluctuations in the condensate, and the faster fluctuations of fi at constant No caused by number-fluctuations 
in the excitations. 

The importance of number-fluctuations in the condensate, assumed at first in the phenomcnological approach due 
to the importance of No for the value of the chemical potential, but later born out by the microscopic calculations, 
leads to the appearance of the linear relaxation-rates 7 C of the condensate-number as an important characteristic 
inverse time-scale of the problem. At times much shorter than j~ x phase-diffusion of the condensate-phase due to the 
fast number-fluctuations in the excitations can occur. In the same regime may also occur the process of collapse due 
to the reversible spreading of the phase caused by the static uncertainty in N and the associated chemical potential. 

At times large compared to 7" 1 the number-fluctuations in the condensate are dynamical and irreversible, and lead 
to the replacement of the reversible collapse by an irreversible phase-diffusion with a larger diffusion-rate than in the 
short-time regime. 

The second and larger part of this paper was devoted to microscopic theory. First we have provided a microscopic 
derivation of the phenomcnological Langevin equation, established microscopic formulas for all phenomenological 
parameters and also exhibited the relation between the short-time diffusion rate and fluctuation rates of the population 
numbers of excitations via a sum-rule. Then the microscopic theory was used to evaluate the transport-parameters and 
the various rates as a function of temperature, particle-number and the scattering length of the interaction potential. 
The evaluation was done for two simple cases - the cubic box-like trap, where the form of the condensate mode does 
not depend on N and the thermal cloud penetrates the condensate homogeneously, and the isotropic harmonic trap, 
where the form of the condensate-mode changes with No and the thermal cloud is located preferentially near the 
boundaries of the condensate. The physically important results for both kinds of traps are similar, even though they 
have to differ, obviously, in the details of the scalings with the atom-numbers and the scattering length. 

The calculation of the transport-parameters reveals some interesting physical results: 

- The fluctuations driving the absolute value |ao| and the phase <f> of ao are quite different in strength, those driving 
|ao| being the much stronger ones. The reason for this is a pronounced squeezing of the bath of thermal excitations 
with respect to the instantaneous phase of the condensate. This squeezing reaches nearly 100% for the lowest-lying 
modes, which is the reason that fluctuations of <j> are practically not driven by such modes. On the other hand, the 
contribution of the high- lying modes to the fluctuating forces driving |ao| and <fr is nearly the same (after the obvious 
normalization with |cto|), i.e. there is no squeezing in this (much weaker) contribution to the noise. 

- The cross-correlation between the fluctuations driving |a | and 4> are found to vanish exactly in a real condensate, 
where both the Gross-Pitaevskii equation and the Bogoliubov-Fetter equations are real and all solutions can (but 
need not) be taken real. This can also be understood as a general consequence of time-reversal symmetry: <f> is a 
velocity potential and therefore odd under time-reversal while |«o| is even under time-reversal. Their fluctuating forces 
therefore transform oppositely. In a time-reversal symmetric condensate (no vortices) the cross-correlation between 
an even and an odd quantity under time-reversal must vanish. 

It turns out that the relaxation rate 7 C of the atom-number in the condensate is the largest of the calculated 
rates. In particular it is larger than the collapse-rate and the phase-diffusion rate which, like j c , are proportional 
to temperature in the regime ksT > [i. It is also larger than the decay rates of the lowest-lying collective modes 
70, which might look surprising because at the same time the theory tells us that "f c is dominated by the particle 
transfer-rates between the condensate and the low-lying modes. However, it is clear that j c ought to be larger than 
70 because the condensate couples to all low-lying modes in parallel which increases the number of decay channels by 
a factor proportional to the ratio of the chemical potential and the lowest-lying mode frequency. 

The next largest rate we find is the thermal phase-collapse rate jcoiiapse- It turns out to have the same func- 
tional dependence on T,a,(N ) and N as the decay-rate of the lowest-lying collective modes. I cannot see any 
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fundamental reason for this coincidence and have to count it just as that. Physically the smallness of jcoiiapse/lc 
means that the phase-collapse will not be observable at finite temperature because it can only lead to a decay-factor 
exp(— \{"icoiiapseHc) 2 ) very close to 1 before phase-diffusion takes over. 

Finally, the phase-diffusion rate is the smallest of the rates calculated here. We find the simple nice result that 
the ratios 7 c /7o and 7o/-D0 are of about equal order of magnitude, given by the ratio of (i to the smallest excitation 
energy, which is Hujo for the harmonic and v/ n/m(2irh/ L) for the box-like trap. Instead of 70 we may also take j c in 
these ratios with the same conclusion. like the rate 7 C is found to be dominated by the atom-number exchange 
between the condensate and the low-lying modes. 

This observation actually explains the coincidence of the two ratios we have just indicated and turns them into a 



precise relation: In (3.19) for we put r = which is exact for real condensate-modes and neglect Tq + V , w hich 



comes from high-lying excitations. Then multiplying the resulting expression for D^, with 7 C = t c 1 from (3.17), we 
readily find 



1 

-^D^ ■ 7c = ~i collapse (9.1) 



with "(collapse from ( |3.23| ) again with T' = 0. 



Let us now compare our results with related ones found in the literature. Most closely related to the present work 
in goal and scope is a paper by Jaksch et al p8| on the intensity and amplitude fluctuations of a Bose-Einstein 
condensate at finite temperature, which builds on extensive earlier work by Gardiner and Zoller with collaborators 
(cf. the references given in |28|]). Unlike the present paper it also takes into account trap losses. The theory presented 
in j?8) is based on a conceptual division of the Bose gas into two energy regions called the condensate band and the 
noncondensate band. In this construction the boundary between the two regions is chosen in such a way that the 
noncondensate band is not significantly affected by the mean field of the condensate, while the influence of excitations 
in the condensate band is neglected. Thus the main physical difference of |2^] to the present work is that it neglects 
fluctuations of particles from the condensate mode to quasiparticle modes as well as to very low-lying one-particle 
excitations. 

By contrast in the present work we avoid the division of the energy region into two parts. We find, as we have 
discussed, that the exchange of particles between the condensate and the low-lying modes makes not only an important 
but in fact the dominant contribution to the relaxation rate of the condensate-number and the phase-diffusion rate, 
determining their dependence on temperature, atom number, population of the condensate and scattering length. 

The importance of the particle-exchange between the low-lying excitations for the phase-diffusion of the condensate 
and the number-relaxation rate 7 C was first pointed out in ||, while for the intensity of the number-fluctuations 
in the condensate this had already been shown in Jl4j]. The theory put forward in |6| already proceeded along 
essentially the same lines we follow here, but it had some short-comings which we overcome and correct in the present 
work: The squeezing of the noise from the thermal cloud with respect to the phase of the condensate was briefly 
remarked upon in fel , but was not taken into account in the calculation of the transport coefficients presented there 



and in the formula for phase-diffusion. Moreover, in the conservative part of the Langevin-equation (3.6) Aq/j, was 
replaced by <9Ai ;l (|a:)| 2 )/<9|a:o| 2 in 0, which, on scrutiny, appears questionable when used in conjunction with the 



fluctuation formula (2.27, f2~2g| ) for (ANq). After all, neither Ap/i nor AF are equilibrium quantities. The use of the 



aforementioned relation between them is therefore avoided here. 

Even though in the present paper I have opted for the use of the fluctuation- formulas ( 2.27 , 2.28 ), which in my 
opinion have a firm basis, it is only fair to mention that they are still under debate in the current literature, see e.g. 
jlql . In another recent paper with some bearing on this topic Bergeman et al. |2^| use as equilibrium distribution 
for the condensate number P(Nq) ~ exp[((/i)Ao — jj(15Noa / 'do) 2 / 5 No) / 'ksT], (cf. the discussion after their eq.(21)), 
which implies (ANq) ~ T(Nq) 3 ^ 5 , a result which is rather different, both in the temperature-dependence and in the 



scaling with the particle-number, from the result (S.19) on which our present calculations have been based. It is 
clear that not the method but the details of our results on the dynamics of the fluctuations of the condensate would 
change, if the results on the statics would be changed. Needless to say that a resolution of the theoretical debate 
concerning the correct approach to the statics seems urgent and would be highly wellcome. Vice versa experimental 
results on the dynamics (i.e. on j c and D^) would also help to decide, by applying the theory presented here, which 
of the approaches to the statics of the number-fluctuations in the condensate proposed in the literature describes the 
physics correctly. 

A quantum kinetic theory of trapped atomic gases has also been formulated by Stoof (3^]. In Q the general 
coupled Fokker-Planck equations of the condensate and the excited modes are presented and applied to the kinetics of 
the formation of a condensate. This problem has also been studied by Gardiner and coworkers |3l| as well as Kagan 
and Svistunov |32|, where also earlier work by further authors is quoted. 
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By contrast the present work has focussed on the fluctuations around the equilibrium state of the condensate, after 
it has been formed. However, the application of our approach to the kinetics of the formation of a condensate would 
be an interesting goal for future work. 

Experimentally the rates 7 C and we have calculated should be measurable. The rate 7 C may be observable as 
the relaxation rate of the condensate back to its equilibrium state after creating a non-equilibrium state by a sudden 
small change of temperature via evaporative cooling. The sum of the phase-diffusion rates of two condensates could be 
measured by monitoring the phase-difference between them after it was initially fixed by measurement or preparation 
at a reference-time t = 0. Methods for measuring phase-differences in Bose-Einstein condensates have recently been 
demonstrated J3|-§| ■ It is to be hoped therefore that the phase-diffusion in Bose-Einstein condensates - a fundamental 
process intimately linked to the spontaneously broken gauge symmetry in a finite system - will be measured in the 
near future. 
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APPENDIX A: 



Here we wish to derive the expression ( p. 13 ) foriJo- Using the Gross-Pitaevskii equation (2.4) we put ( |2.8| ) in the 
form 

Ho = (mo - <M»N| 2 - ^Nl 4 f d 3 x^ Q (Al) 



Taking the derivative with respect to | ckq | 2 we get 



9H ° Mo - <M> + Kl 2 ( -£% -U f d*x^ - \a \ 2 U [ ) (A2) 



d\a \ 2 \o\a \ 2 J J o\a 

To evaluate this further we use eq.([2.4|) and its derivative with respect to |ao| 2 



k2 V 2 +V mo + 3C/oKlV 2 ") 4% = ( #^ " Uort ) iPo (A3) 



2m v ™ Ul U| r V <9|a | 2 V^ol 2 



Multiplying eq.(A3) with ipo and integrating over space, using the Gross-Pitaevskii equation (2.4) after partial inte- 



gration, we derive the identity 



'9|« | 2 £>KI 2 

which is used in in (A2) to yield dHo/d\ao\ 2 — Mo — (m)j and upon integration results in ( p. 13 ) 



[/oKI 2 / d\^4% = -U [ ^d 3 x (A4) 



APPENDIX B: 



Here we wish to derive eq.( 4.14 ). This is achieved if we succeed to show that the coupling of the condensate and 
the thermal cloud via 



£3 = U ^(Nj) J d 3 xM + (e- 1 ^ + e^"x)x (Bl) 
gives rise to the systematic change of 3(£(i)), to first order in the interaction, of 



Wt®)* = - 2 4Mf dt ' S Mt~t')^, (B2) 
nk B T J_ 00 o\a \ 2 
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because this can then be used in (4.7) to yield (4.14). In (Bl) we could put |cko j — y (^o) since we linearize around 
equilibrium and only wish to calculate the dissipation in |oop which is conjugate to <fi, the variable we kept in (Bl). 
Standard first order perturbation theory with adiabatic switch-on of the interaction gives, with e — > +0, 



We can rewrite this as 

Smi(t)))^ = ~2iVm f dt' (^{trfe-W +x'; f+ (M'K 0(t,) ) e et ' 

J — oo 

where we introduced the response functions 



(B3) 



(B4) 



with 



y!' Jt (t, 0= <[»(&))> £ + (*')]>*©(* - 



|(t) = c/ / d 3 ^* (*)x(*)x(*)- 



(B5) 



(B6) 



Here Q(t — t') is the Heaviside step-function. We shall define 0(0) = without loss of generality. The fluctuation- 
dissipation theorem (in the classical frequency domain hu> <C fcsT) ensures the relations 



2fc B T dt' J « v 



" ,/ i@(t-t') d , 



with the correlation functions 



2k R T dt' J « + 



s^(t,0=<3(e(*))l(* , )>* 



s Ji+ (t,t')={Z(£(t))t + (t'))t. 



We can use (|B7j) in (B4) and apply a partial integration in t' to obtain 



= f dt' (Sj^t, t')e-^ - S Ji+ (t, t')< 

1KB 1 J —oo 

(^(M)e-^ +^- + (M)e^) 



#(0 



which can be rewritten as 



dfS. 



• 7j(M) ^ 



(B7) 



(B8) 



(B9) 



(BIO) 



11 9(^(t)) according to eq.(4.9) contains an explicit external time-dependence via <fi(t), in addition to the internal time- 
dependence of x(t), x + (t) via their Heisenberg equations of motion. This explicit time-dependence has to be taken into account 
when applying the fluctuation-dissipation theorem. We avoid this additional step by applying the time-derivative in the 
fluctuation-dissipation relation (B7) directly to the second time argument t' , of course with the appropriate extra minus-sign. 
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with Sjj(t,t') defined by eqs.(gj. 



, 1111 ) and 



(f,0 = ^((|(t)e _i * (t) -| + (t)e^W)(|(t')e-^ (t,) +^ + (0e i * (f ' ) )) f 



(Bll) 



The constant term with Sjn(t,t) amounts to a small shift of the equilib rium value of fi in the final result which we 
shall negle ct like other terms contributing to suc h shif ts. Then using eq.(4.2) we put Ud(j){t')/dt' = — dHo(t')/d\ao\ 2 
in eq.(B10|) which establishes (Bl) and hence eq.(4.14). 
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